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Abstract 

When studying the production of minijets and transverse energy flow in 
high energy hadron-hadron or nucleus-nucleus collisions, two essential points 
have to be taken into account. First, one has to account for the virtuality of 
the colliding partons and secondly, it is important to avoid double counting, 
when many links in a parton chain can be interpreted as the momentum 
transfer in a hard subcollision. The Linked Dipole Chain model, introduced 
for low-x DIS, is particularly suitable for a study of these problems. It 
describes (mini)jet production in a /^-factorizing formalism, which includes 
all links in a parton chain on an equal footing, avoiding double counting. 
In a "naive" calculation based on integrated structure functions, the cross 
section blows up for small p±, which makes it necessary to introduce a soft 
cutoff. In our approach we find a dynamical suppression at low p±, which 
makes it possible to extrapolate to higher energies and make more reliable 
predictions for RHIC and LHC. 
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1 Introduction 



With increasing energies in hadron-hadron or nucleus-nucleus collisions, the cross section 
for hard subcollisions increases and becomes more and more important. The amount 
of minijets and transverse energy becomes essential for understanding the background 
in searches for new particles or new phenomena at the LHC as well as for the "initial 
conditions" in nucleus-nucleus collisions. In calculations of the flow in either a quark- 
gluon plasma or a hadronic phase, the results are very sensitive to the properties of the 
initial parton state which resulted from a large number of hard subcollisions. Thus a 
reliable estimate of these initial conditions is essential for the interpretation of signals 
from a possible plasma formation in experiments at RHIC or LHC. 

At high q± the jet cross section can be described by a product of structure functions 
describing the flux of partons and the cross section for a hard partonic subcollision. 
Symbolically we write (cf. Fig. |l]a) 



This is a relevant description when q± is so large that the structure functions can be 
described by DGLAP evolution, i. e. by /cj_-ordered chains from the incoming hadrons 
towards the hard subcollision. In the hard collision cross section, da/dq\, the colliding 
partons are then approximately on shell. For smaller q\ and/or larger energies (which 
implies smaller x-values) we enter the BFKL region, in which non-ordered chains are 
important. This has two essential consequences: 

i) The virtuality of one or both of the colliding partons {k 2 ±1 and/or k\ 2 in Fig. [I]a) can 
be larger than the momentum transfer q\. This implies that the virtualities of the colliding 
partons have to be taken into account. This can be achieved by means of nonintegrated 
structure functions and off-shell subcollision cross sections. 

ii) In a chain with non-ordered fcj_-values a single hard subcollision cannot be isolated. 
There may be several links in a single chain, which can be regarded as a hard collision. 
This situation can be analyzed in a formalism in which every link in Fig. [I] can be re- 
garded as a hard subcollision. This implies that each final state parton is active in two 
different subcollisions. (The number of links in the chain determines both the number 
of subcollisions and the number of final state partons.) Therefore special care has to be 
taken to avoid double counting. 

A formalism for DIS which interpolates between the DGLAP and BFKL regions, and 
which describes both the cross section and the final state properties, was developed by 
Ciafaloni, Catani, Fiorani, and Marchesini fl|]. This or equivalent formalisms can be used 
e. g. to calculate the production of heavy quarks in hadronic collisions. In this reaction 
there is a single link which contains the heavy quark, and the CCFM formalism or the 
semihard approach can be used to calculate the unordered evolution towards the heavy 
quark from the projectile and the target ends f2[. Thus these formalisms can be used to 
solve problem i) above. 

For non-Zc^-ordered chains the single diagram in Fig. p] can correspond to many hard 
subcollisions, e. g. 1 + 2 — > 3 + 4 and 5 + 6^7 + 8 in Fig. [I|b. For the first reaction, 
1 + 2 — > 3 + 4, the section of the ladder between partons 2 and 5 should be regarded as an 
evolution from the top (projectile) end, while for the reaction 5+6 — > 7+8 the same section 
should be regarded as an evolution from the lower (target) end of the chain. A description 
of this process is most simple in a formalism in which the evolution is explicitly symmetric 
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Figure 1: A fan diagram for a hadron-hadron collision. 

between the two ends. It is also essential that the formalism can be interpreted in terms of 
production probabilities for exclusive final states. The symmetry is not a trivial feature 
for the following reason. The fan diagrams in Fig. [I] should be regarded as the initial 
state radiation. To get the complete final state, final state emission has to be added in 
specified kinematic regions. The separation between initial state and final state emission 
is not determined by Nature, but is defined by the calculation scheme, and therefore the 
regions for final state emission depend on the formalism used. In the CCFM approach this 
separation is not symmetric, and consequently also the initial state ladder is asymmetric. 
Thus the CCFM formalism is not immediately convenient to describe many subcollisions 
in a single chain, and thus to solve problem ii) above. 

The final state emission, when one gluon is split into two, should not significantly 
affect the flow of energy. Note, however, that the distribution of produced gluons is not 
a uniquely defined quantity. Besides the arbitrariness in the separation between initial 
and final state emission, discussed above, also the final state radiation depends on the 
resolution chosen in the definition of an individual gluon, i. e. on the cut-off in the final 
state cascades. 

The Linked Dipole Chain (LDC) model |J is a reformulation and generalization of 
the CCFM formalism. The separation between initial and final state radiation is chosen in 
such a way that the description is explicitly symmetric between the projectile and target 
ends of the ladder. It can within a single formalism describe different types of events, 
"normal DIS", boson-gluon fusion events, and events with hard subcollisions in resolved 
photon events. For such a resolved photon event the result can be directly interpreted 
as evolutions from both the proton and the photon ends towards a hard parton-parton 
collision. This feature implies that the LDC model is also particularly suitable for a 
description of hard collisions in hadronic interactions. For any link in the ladder, the 
contribution can be directly interpreted as evolution from the two ends multiplied by an 
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off-shell parton-parton cross section for the link under study. This symmetry property of 
the LDC model makes it useful in solving problem ii) above (the double-counting is easily 
corrected for in this formalism). 

It is well known that the BFKL formalism with a constant coupling, a s , implies 
that the transverse momenta grow like a random walk in y or in In 1/x (what is called 
Bartels' cigar). For a running a s we obtain instead a saturation of the fcj_-distribution 
Q. Therefore very large chains develop a central plateau in rapidity. This means that 
with increasing beam energy the transverse energy density in each chain will stay limited, 
although the total dE^/dy will increase as the number of possible chains grows with 
energy. The fact that the /c^-distribution stays limited also for high energies implies that 
the result is sensitive to soft nonperturbative physics, and thus depends on a necessary 
cut-off for small k±. This cut-off can, however, be fixed by a phenomenological fit to 
experimental data for ^(x, Q 2 ), and therefore does not imply a large extra uncertainty. 

When comparing our result with the "naive" expression in eq. (|I|) we find that the 
latter significantly overestimates the cross section for smaller q±. In order to agree with 
experimental data for E±-Qow, many calculations based on eq. (|l|) introduce a cutoff for 
small q± of the order 2 GeV/c 0. This implies, however, that it is difficult to make 
predictions for higher energies. Without a dynamical understanding of the origin of the 
overestimate, it is not possible to judge how an effective cut-off should vary with energy. 
The dynamical suppression in our formalism corresponds to such an effective cut-off, 
indeed of the order 2 GeV/c, which is growing slowly with the total collision energy. 

The results in this paper are mainly based on approximate analytic calculations, which 
demonstrate the qualitative features. In a future publication we want to present more 
detailed results based on MC simulations, where e. g. corrections from quark lines, 
subleading terms in the splitting functions and exact energy-momentum conservation, are 
also taken into account. 



2 DIS 

A deep inelastic scattering event is generally described in terms of a fan diagram as shown 
in Fig |2|. Here denote quasireal partons emitted as initial state radiation, while the 
links hi are virtual. The dashed lines denote final state emission, which is assumed to be 
emitted without changing the cross section and with negligible recoils for the emitting 
partons g^. 

In the large Q 2 region, the DGLAP region, the dominant contributions come from 
ordered chains which satisfy Q 2 > k\ n > ... > k\^ > k\^_ x and k +)i > Each such 

chain gives a contribution (xj = k +)i /P +)tot ) 

J\ _ dxi dkj_ j _ 3a s . , 

\\ol -y^- where a = . (2) 

Integrating over the appropriate integration regions and summing over possible values of 
n, the number of links in the chain, we readily obtain (for a fixed coupling a) 

F(x,Q 2 ) ~]>>" (lnl ( x)n(ln ^ 2)W « exp(2 v /«lnQ 2 lnl/x). (3) 

For a running coupling we get instead of InQ 2 a factor lnlnQ 2 . We want to stress that 
to get the properties of the final state, we have to add final state emission within regions 
allowed by angular ordering. 
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Figure 2: A fan diagram for a DIS event. The quasireal partons from the initial state 
radiation are denoted qi, and the virtual propagators k{. The dashed lines denote final 
state radiation. 

For very small x and limited Q 2 , the BFKL region, also non-ordered chains become 
important, although suppressed. Solutions to the BFKL equation || increase like a power 
l/x x for small x- values. Such a power-like behavior with A ~ 0.3, is indicated in data on 
F 2 from HERA. (We note, however, that it is also possible to describe this increase by 
NLO DGLAP evolution.) 

In the interpolation region between the DGLAP and BFKL regimes we have to cal- 
culate suppressed contributions from non-ordered chains. For each chain of initial state 
radiation, final state emission should be added within specified kinematic regions. This 
final state radiation should give negligible recoils. It should also be described by Sudakov 
form factors, which implies that the final state emission only affects the properties of the 
final state and not the cross section (i. e. the reaction probability), which is described 
by the structure function. Thus a chain with specified initial state radiation represents 
a set of final states with all possible final state emissions within the allowed regions. We 
want to stress that the separation between initial and final state radiation is not given by 
Nature, but is defined by the calculation scheme. 

A specific scheme, which interpolates between the DGLAP and the BFKL results, 
was presented by Ciafaloni, Catani, Fiorani and Marchesini, the CCFM model Q. In this 
scheme those final state gluons, which are not followed in rapidity (or angle) by a more 
energetic gluon, are regarded as initial state radiation, all other as final state emission. 
With this definition they showed that for small x the nonintegrated structure function, 
/, is given by the expression (see ref U for details) 

f(x, kl, q 2 ) ~ £ / a — ^A ne (z t , k 2 ±ii , q 2 )5(x ~ n Zi )S(kl - kl n )5(q 2 - q 2 ). (4) 

Here q± t i and k± i are the transverse momentum of the real and virtual partons as indicated 
in Fig. ||, and the variable q is defined by the relation q = q±/(l — z), and the function 
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A ne (z, kj_, q 2 ) is a specific non-eikonal form factor. The integration region is restricted by 
the relation 

q±,i+i > Ziq± ;i /(1 - Zi) (5) 

which follows from angular ordering, and by the kinematic constraint (called the consis- 
tency constraint) 

k ±,i > z iQ±,i- (6) 
We note that this nonintegrated structure function f(x, kj_, q 2 ) depends on two scales, k± 
which defines the transverse momentum of the last link in the chain, and q which depends 
on the angle of the last emission, and therefore specifies the boundary of the angular 
region in which later emissions are allowed. 



3 The Linked Dipole Chain Model 

The Linked Dipole Chain model || (LDC) is a reformulation and generalization of the 
CCFM model. In the LDC model more gluons are treated as final state radiation. The 
remaining (initial state) gluons are ordered both in q + and in q_. (This implies that they 
are also ordered in angle or rapidity y.) Thus a single chain in the LDC model represents 
a set of chains in the CCFM scheme, all with the same "backbone" of harder gluons. In 
ref it is demonstrated that if we sum over all states in this set, with their corresponding 
non-eikonal form factors, then all this adds up to unity. Thus the form factors are exactly 
canceled, and a nonintegrated structure function T can be written in the simple form 

JJa — - — ^9(q +)i -i - q +t i)8(q^ ti - - Uzi)5(\nk^ - lnk 2 ±n ). 

(7) 

We note in particular that this result is symmetric in q + and g_. (An essential point 
for this symmetry and the cancelation of the non-eikonal form factors is the consistency 
constraint in eq. @.) The fact that the ordering in both q + and g_ also automatically 
implies an ordering in y (i. e. in angle), means that the nonintegrated structure function 
J-'lx, k'jj depends only on one scale, k\. (In contrast, the corresponding form factor, 
f(x,kj_,q 2 ), in the CCFM model also depends on the scale q 2 , related to the boundary 
of the angular region allowed for further emissions.) This fact implies a considerable 
simplification of the formalism. The symmetry also implies that the chain in Fig |] can 
be interpreted either as evolving from bottom to top (i. e. from the proton end towards 
the photon) or evolving from top to bottom (from the photon towards the proton). This 
feature will be essential for our analysis of the minijet distribution. 

Fig. ^| shows a typical chain in a (y, lng^)-plane. The real emitted partons qi are 
mapped onto points in this figure. The lightcone components In q± = ±y + 1 In q\ grow in 
the upper right and upper left directions as indicated in the figure. The virtual propagators 
do not have a well defined rapidity; they are represented by horizontal lines, whose left 
ends correspond to the values of k±^ and k +> i, while the right ends correspond to k±^ and 
fc-^. The larger region allowed for final state radiation in the LDC model corresponds to 
the region below the horizontal lines in Fig. [3]. 

It is convenient to express the result in terms of the propagator momenta k iy using the 
relations d 2 q± ti = d 2 k ±!i and q\ { ~ max(/c^ i , k\ Thus eq. (|7|) can be written in the 
form 

•^~E/n«— — T^hi y (8) 

J Zi max^fcj_ j, j i J 
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where we have suppressed the 9- and ^-functions. This result implies that for a "step up" 
or a "step down" in k± we find the following weights 



d 2 q±,i d 2 k ±ti 

2 ' ~ -T2 — > k±,i > fc_L,i_i and (9) 

— 2 ~ -T2 , 2 «L,i < «±,i-l- (10) 

<7±,i fc _L,i fc _L,i-l 

Thus for a step down we have an extra suppression factor k 2 Li /k 2 Li _ l . This implies that 
if the chain goes up to k± )max and then down to k±j ina i we obtain the factor 

file 2 k 2 

n uhv ±,i final 

~^'~k^~- (U) 

" / _L,i rh ±,max 

This gives a factor l/^x mm ., which corresponds to a hard parton-parton subcollision. 

In DIS also chains for which k\ j inal > Q 2 can contribute to the cross section. These 
chains correspond to boson-gluon fusion events, and the definition of the structure function 
F(x, Q 2 ) implies that these contributions contain a similar suppression factor Q 2 / k\ j inal . 
Thus the relation between the integrated and the nonintegrated structure functions, some- 
times symbolically written as F(x,Q 2 ) ~ / ^-J-'ix, kV), is more explicitly given by the 
relation 

F(x, Q 2 ) = f ^{x, k 2 ± ) + J q2 d -j^F(*§, ■ (12) 

In the second term, besides the suppression factor Q 2 jk\ we have also a shifted x- value. ' 
'Normal" chains, for which k\ ^ inal < Q 2 , end on the line AB in Fig. f|, which corresponds 
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ln(Q 2 ) ln(l/x) 



Figure 4: Different types of parton chains. I: "Normal" DIS. II: Boson- gluon fusion. Ill: 
Hard resolved photon collision. 

to k + j ina i = x ■ P+,tot, but due to energy-momentum conservation, chains for which 

k 2 

k±jinai > Q 2 en d instead on the line AC, which corresponds to k + ji na i = x ■ ± ^ 2 nal ■ P +)tot . 
Note that the relation in eq. fll~2| ) is a leading approximation valid when k\ is either much 
smaller or much larger than Q 2 . Exact energy- momentum conservation implies that the 
argument in T should be x(k"j_ + Q 2 )/Q 2 , which approaches x or xk\jQ 2 in the two limits. 
The exact kinematic relation is included in the MC in the whole chain including this final 
vertex. Our experience indicates that the approximations presented here are sufficient to 
understand the essential properties. 

We see that in the LDC model different types of reactions are treated in the same 
formalism, without double counting or missed parts of phase space. Thus Fig. |] shows 
three chains representing different types of reactions: 

I. "normal" DIS, 

II. boson-gluon fusion: k\ t incd > Q 2 , 

III. hard resolved photon collision: k\ max > k\j inal {> Q 2 ). 

As mentioned above, the LDC formalism is fully left-right symmetric, meaning that 
the same result is obtained if the chain is generated from the photon end instead of from 
the proton end. Although not evident from eq. (|S|), this is obvious from the expression 
in eq. (0). This feature also means that the same formalism can be used if the (resolved) 
photon in one end of the chain is replaced by a hadron. In section 5 we will see that the 
symmetric property of the LDC formalism makes it particularly effective for an under- 
standing of the minijet distribution in hadronic reactions or nucleus-nucleus collisions. It 
is also important that the result in eq. (0) not only describes inclusive properties, but can 
be interpreted as the production probability for an exclusive final state. 

Before we go into more details about the properties of this minijet plateau, we want 
to mention some other properties of the LDC model: 

• It is straight forward to allow for a running coupling a s within the formalism. 
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• It is possible to include quarks and other non-leading effects (non-leading terms in 
the splitting functions and exact energy-momentum conservation). The result is also 
improved when the link with highest p± is adjusted to the exact matrix element. 

• The chains contain only those gluons (or quarks), which cannot be regarded as final 
state radiation. This might be called the "backbone" of the chain, and the fact that 
hard jets are not yet subdivided into many sub-jets may make it more easy to study 
minijets and E^-flow, and to interpret the results of the calculations. 

• The fact that there are fewer gluons in the primary chain also implies that typical 
^-values are smaller, and therefore smaller sub-leading effects are expected. 

• The formalism is suitable for MC simulation. Such a program is developed by 
Lonnblad and Kharraziha 0. 

Naturally it is essential to verify that the LDC model also can reproduce experimental 
data. Preliminary results indicate that the LDC MC indeed is able to successfully describe 
experimental results, both for the structure functions and for the properties of the final 
states, for example the production of jets and transverse energy flow 0. 



4 Integral equations and asymptotic behavior 

It is straight forward to derive integral equations for the non-integrated structure function 
T . From the relation in eq. ( |T2D we obtain the equation 

df(!, re) 



dl 



[ K dK'a(K)F(l, re') + f d^a{^)T(l + re - re', Kf)exp[-(*f - re)] (13) 



where / = ln(l/x) and re = \nk\. 

It is well known that the BFKL formalism with a constant coupling, a s , implies 
that the transverse momenta grow like a random walk in lnk±, which gives a Gaussian 
distribution that widens with ln(l/x), <lnAr[>~ \J\n.{l/x). A running coupling favors 
smaller /c_i_-values, and this implies that the /^-distribution does not widen indefinitely, 
but saturates for small x |3| . The solution to eq. flllf ) can then be written in the factorized 
form 

T ps — • /(re), x small. (14) 

The fact that small k±- values are not suppressed implies that the evolution and the value 
of A is sensitive to the soft region, and therefore cannot be determined from perturbative 
QCD alone. (This feature is consistent with the very large lower order corrections to 
the BFKL equation 0, which indicates that the perturbative series may converge badly.) 
This implies that some kind of cutoff, Qo, has to be introduced, and determined from fits 
to experimental data. If this cutoff is determined from fits to data on F 2 from DIS, then 
this result can be used in our calculation of the minijet distribution in hadronic collisions. 
Although the dependence on x and k± factorizes as expressed in eq. flI3|), the two factors 
are not independent. The relation in eq. ( |I~3"D implies a correlation such that the power A 
also determines the /c^-dependence. Thus for large k±, f takes the asymptotic form 

/(re)ocre^- 1 (15) 



8 



where a is defined by the relation 



3a s _ a 36 

a = — = HQW) ao = 3Y^-/ (16) 



Actually MC calculations show that the simple power in eq. fllq) is a surprisingly good 
approximation, not only for large values of k± (or k = \n(k"j_/A 2 )) but for the whole 
fc^-interval [§]]. 

The factorized form in eq. QT3D implies also, that for limited Q 2 and very small x (or 
hadronic collisions at very high energies) a central plateau is developed in the minijet 
distribution. The properties of this plateau will be further studied in the next section. 



5 Inclusive Jet Cross Section 

In this section we present asymptotic results and general properties. A more detailed 
analysis will be presented in a future publication. 

The properties of the fan diagram in Fig ^| depend on the chosen separation between 
initial state and final state radiation. As mentioned above, the same event corresponds to 
a chain with fewer links in the LDC formalism than in e. g. the CCFM formalism, because 
in the LDC model more emissions are formally treated as final state radiation. The number 
of jets is not a well-defined concept, if it is not accompanied by some specification of the 
resolution, and two different schemes can be equally correct if one jet in one scheme 
corresponds to two or more sub-jets in the other. The total -E^-flow is however expected 
to be approximately unaffected, when one jet is split in two or several smaller jets by final 
state radiation. As mentioned in section 3, the LDC formalism is fully symmetric with 
respect to the two ends of the fan or ladder diagram, and we will see that this makes it 
particularly effective for an analysis of the minijet distribution. (In the CCFM formalism 
the fan diagram for the initial state radiation is not symmetric; the symmetry is restored 
only after inclusion of the final state emission.) It is also important that the formalism 
can be interpreted in terms of production probabilities for exclusive final states. 

A long chain as in Fig. |5|, with a soft probe (Q 2 ~ Qq) and with a local maximum 
k±, m ax = k± t i, corresponds to a hard scattering between two partons with momenta 
and —k i+ i. Due to the symmetry of the expression in eq. (|7|), these two pieces of the 
chain can be interpreted as evolution from both ends towards the central hard scattering. 
The momentum transfer in this hard sub collision is given by fcj, and there is a minus sign 
in the momentum of one of the colliding partons, as this part of the chain is regarded as 
evolving from the top of the fan diagram, downwards towards —k%+i. 

A single chain may also have more than one local maximum, which corresponds to 
two or more correlated hard sub collisions. If we want to obtain the inclusive jet cross 
section, or the total i?i-flow, we have to include all produced partons qj in the chain. To 
calculate this we start by studying a single link in the central part of a long chain. There 
are three different types of links. Besides a local maximum as in Fig. |5| and Fig. ^|a, we 
have also the two possibilities shown in Fig. |6|b and Fig. ^|c. To simplify the notation we 
call the gluons k a , k, and fc& as indicated in Fig. ^, and study the three cases: 

1. k± > k± ta ,k_L,b- Transverse momentum conservation implies that q±, a ~ qj_,b ~ k±, 
and it corresponds to a normal hard Rutherford scattering between two partons. 
The two colliding partons have virtualities given by k\ a and k\ b , which both are 



9 




Figure 5: A local maximum transverse momentum, k± )max = k± t i, in a long chain, corre- 
sponds to a hard subcollision between partons with momenta ki-x and — fcf+i. Xi and x\ ^ 
are defined as x^ = k + ^/P + ^ ot and x\~^ = —k_ i/P_j ot respectively. 



small compared to the exchanged momentum k\. From eqs. (|3|-|l0|) we see that the 
link is associated with the following weight factor 



1 

k 



, <(ki), (17) 



which corresponds to the cross section for gluon exchange between two quasireal 
particles. 

2. k± t b > k± > k± >a , which implies that q±^ a ~ k±_ and q±j ~ k± t b- From eq. (|S|) we see 
that instead of the expression in eq. (|17D , this link is associated with the following 
weight 

' a s (kl)-a s (kl b ). (18) 



h 2 ■ h 2 



A similar result is obtained for k± }Q > kj_ > k± t b- 

3. k± < k± >a , in which case g_L jQ ~ k±^ a and q± t b ~ k± t b. The corresponding weight 
is now given by 

' ■<kl.)'a B (k 2 , b ). (19) 



h 2 ■ h 2 

K ±,a li ±,b 



We see that when k± is less than kj_, a or &_l,6> one factor l/k\ is replaced by l/k 2 ia 
or l/k 2 Lb respectively. As will be discussed more in the following, this implies that the 
inclusive cross section becomes non-singular for small transverse momenta. (As discussed 
after eq. (|12"D non-asymptotic contributions to these weights can be included most easily 
using the MC.) 

Keeping k a , k, and k b fixed, we note that integrating and summing over all par- 
tons to the right of the link corresponds exactly to the non-integrated structure function 
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Qb 



(b) 




Figure 6: The three different possibilities for a link. 



J-(x a , k\ a ), where x a is the Bjorken variable, i. e. the scaled positive lightcone momen- 
tum, x a = k +ta /P +ttot . Similarly, integration and summation of partons to the left of 
the link gives Ffab, k\ b ), with x 6 the corresponding scaled negative lightcone momentum 

X b = x b ^ = ~k-,b/ P-,tot- 

The inclusive jet distribution can then be written in the form (s is the Mandelstam 
variable for the hard subcollision) 



da. 



incl 



dq]_dy 



x 



x n 



dxb 
x b 



dk\ 



k 2 



3~(x a , kj_ a ) ■ J-(xb, kj 



1 fl!L(, - / 2 ■' 2 

2 dq\ 



(q±-> k±^ a , k^ h i s — x a Xbs) 



1 x 
8{y--\n—). 

2 x b 



(20) 



In this formalism each produced parton is counted in two different links, one to the left 
and one to the right. In our notation J^- is the inclusive cross section counting both 

emitted partons, and therefore a factor | is needed to avoid double counting. The cross 
section is obtained by integrating eqs. ( |T7l - |i9l) over k\ with the constraint k\ < s. In 

case 1 above we will then get the integral (for k±^ > k±,a) 



dk 2 ■ ^. ■ 2 • %i - kl) = - ■ 9(ql - k 2 ^) ■ 9{s - ql) 



k'\ 



(21) 



Here the factor 2 in front of the ^-function originates from the fact that we count both q\ a 
and qj_ b , each being approximately equal to k\. The contributions from cases 2 and 3 can 
be calculated in the same way. The result obtained for a running coupling is presented 
in the appendix. The qualitative features can more easily be understood from the result 
obtained for a constant a s , in which case we find (again for k± t b > &±,a) 



da 
dq\ 



oc a. 



+ 



+ 



- " %1 - Kb) ■ 0{s - ql) 



l 



q\ ■ k] 



0(Ql-kl, a )-9(kl b -ql)+5(ql-kl b 



S(ql ~ k") + S(ql - k] 



A i 



max - 
\ s 



■In 



Ql 



Ay I th I 



,a _L,t 



k 2 



+ 



(22) 



Here the first line corresponds to a "normal" hard subcollision, case 1 above. The second 
line corresponds to case 2, where one of the colliding partons has a virtuality larger 
than the exchanged momentum —i=kl, and the third line to case 3 where both initial 
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partons have high virtualities. Note that if we have a soft cutoff Qq, the lower limit in 

k 2 k 2 

the /^-integral f° r case 3 is given by max( ±,a . ±,b , Ql). 

The inclusive jet cross section is now readily obtained from inserting eq. (|22]) into 
eq. (p0|) and integrating with respect to x a ,x b , k± >a and k± tb . Due to the factorized form 
T{x,k) ~ x~ x -f(n) in eq. (|14|), with a power-like dependence on x, the result is inde- 
pendent of the rapidity of the pair, y = |ln|^-, for fixed value of x a -x b = s/s. Thus 
the distribution corresponds to a central plateau with a height proportional to (xi-x 2 )~ x , 
which for fixed s grows with energy proportional to s x . 

The integration is straight forward if we assume the power-like approximation for /(k) 
in eqs. ([H|-[n|), T{x, k) « x~ x ■ f(n) oc x~ x k~$~~ 1 . As an example, for the term in eq. fl2"T| ) 
(the first term in eq. (0)) the integral over x a and x b gives 

^^(y-i]*^).-!/ -{-) X = \(^\ (23) 
x a x b 2 x b {x a x b ) A J q \ s \sj A \q±J 

while the integral over k^ a and k^ b gives (including also the symmetric case k\ a > k\ b ) 

I 2 



dK b ■ n b x / ck a • K a x = a 2 s (q] 



a 



(24) 



The full result for a running coupling, presented in the appendix, is rather lengthy and 
we write it in the form 



The factor s x /q 2X is a consequence of the increase of JF ~ x~ x for small x (cf eq. (|T4])), 
and we have also extracted a factor aKqjJ- Thus the function h(q 2 ±) is defined in such 
a way, that it would be constant in the unrealistic case, where the parton flux is given 
by a scaling function F(x) ~ l/x x , which is independent of Q 2 . The result in eq. (|25|) is 
proportional to s x , but we emphasize that this is the growth of the cross section for jet 
production, and not the number of jets per event. 



Results 

For a quantitative estimate of the jet cross section and a comparison with the "naive" 
approach, it would be suitable to use a Monte Carlo, which can take into account also 
effects from non-asymptotic energy, non-leading terms in the parton-parton cross section 
and contributions from quarks. However, we believe that the qualitative features can be 
understood from an approximate analytic calculation, provided the same approximations 
are used in both approaches. Thus for F(x, qj_) and JF(x, kj_) we use the relations in 
eq. ( |T2"D and eqs. ( jl~4^|l~5|) . We study purely gluonic chains and for the subcollision we use 

the leading expression JS- oc s g 4 X . For the parameters A and Qo (the cut-off for small 
kj_) we use the values 0.3 and 0.6 GeV/c respectively, obtained from the MC fit to HERA 
data. Since we study purely gluonic chains, we may take rif = 0, which implies that the 
power ^ — 1 in eq. ([15]) is close to 2. The result, presented in the appendix, can be 
written in terms of simple elementary functions when this power is an integer. Numerical 
calculations show that the result is insensitive to the exact value of the power, and that 
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the result using ^ — 1^2 agrees well with results from the MC. For these reasons we 
use this value in the results below, and thus assume the following form: 

T{x, k) oc x~ x ■ k 2 (26) 

which implies that the integrated structure function has the form 



F(x, k = InQ 2 ) oc x 



(27) 



(Note that with our convention F(x) describes the density in \nx, and thus corresponds 
to F 2 . Since we have a hard interaction probe we have included the running of a s in the 
final vertex, which implies an extra factor InQ 2 / lnk± in the last term in eq. (12).) 

The transverse momentum dependence of the function h(q±) for a running a s is shown 
by the solid line in Fig. [7[ Also shown in this figure is the contribution from "normal hard 
sub collisions", i. e. partons associated to a local maximum in the k± chain. We see that 
this contribution dominates for large qj_, while for smaller qj_, the partons which are not 
associated to a local maximum play a more important role, contributing more than 40% 
for q\ < 5GeV 2 . For large qj_, h(qjj behaves as a power of lng 2 , which corresponds to 
the scaling violation in the structure functions. For smaller qj_, h(q±) is essentially linear 
in q\. This implies that the total transverse energy flow dE±/dy = ^ ■ / Jp d dq\ ■ q± 

is convergent for small q± (for A < 0.5). Consequently the E± distribution is limited also 
without a low q± cutoff. In Fig. |7|a we also show the g_i_-dependence obtained from the 
MC (normalized to the same total flow), and, as mentioned above, we see a very nice 
agreement with our analytic result. 

We note that although "normal hard collisions" dominate for large q±, the other 
contributions correspond to approximately 25% even for q 2 ± « 1000 GeV 2 . We expect 
that most of this difference would be accounted for in a calculation where hard collisions 
are calculated to next-to-leading order, including also 2 — > 3 parton reactions, which thus 
can take into account one extra parton. This approach would, however not solve the 
problems encountered for small or medium q±. 

In the literature the inclusive jet cross section is often estimated from a product of 
two (integrated) structure functions F(x, g 2 ) and a hard subcollision cross section ^p- for 
quasireal colliding partons. This estimate can be written in the form 



oc ( d J^ d J^5{y-\ln X ^).F{x a ,ql).F{x h ,ql).\ 
dq±dy J x a x b 2 x b - q]_ 



(2c 

<?_L 



Here the function h naive is defined analogously to the function h in eq. (|25|) . 

It is interesting to compare this "naive" estimate with our result in eq. (|25]) . The 
relation between F(x, qj_) and the non-integrated structure function T(x, kj_) is shown in 
eq. fll2"D - There are two contributions to F, one (ordered) contribution with q±_ > k±, 
and another suppressed contribution with q±_ < k±. When one or both colliding partons 
are virtual, there is a suppression factor q\/k\ a and/or q\jk\ h . With a cross section 

2 / 2 \ 

assumed to be proportional to " a 4 as in eq. (|25|) , we see that each link is given just 
the weight presented in eq. fll7D , eq. (|i~8|), or eq. (|i~9P . The difference lies in the fact 
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ql (GeV 2 ) q{ (GeV 2 ) 

(a) (b) 

Figure 7: The transverse momentum distribution of h(qj_) for (a) q\ < 10 GeV 2 and 
(b)q\ < 1000 GeV 2 according to analytical calculations (continuous curves) and MC 
simulation (short-dashed line). (The error bars are the estimated statistical uncertainties.) 
Also shown is the contribution from "normal hard subcollisions" (long-dashed line). The 
scale is arbitrary. 



that in eq. ( p8| ) both outgoing partons are assumed to have transverse momenta given 
by the momentum transfer in the collision. This underestimates the q± for collisions 
corresponding to the links in Fig. |b and Fig. ||c. This underestimate is compensated by 
the fact that every link in the fan diagram corresponds to a hard collision in eq. (|28|) . 
Thus this equation corresponds to two produced partons per link instead of one. This 
double counting does not give a full factor of two, however, due to the underestimate of 
the q± mentioned above. The result is illustrated in Fig. The really emitted partons in 
a chain are marked by a dot, while the naive expression in eq. ( f28|) corresponds to a final 
state parton at every point marked by a circle. Thus in the naive calculation only the 
partons at a local maximum are given their proper weight. Those at a minimum value 
for k± should not be included at all, as no real emitted partons have these fcj_-values, and 
those in-between are counted twice. 



In q\ 



Figure 8: The outgoing partons according to a correct formalism (dots) and according to 
a naive approach (circles). 

Fig. H shows the ratio htq'jj / h naive (q 2 L ) as a function of the transverse momentum 
(note the logarithmic g^-scale). For large q± the Rutherford contribution dominates in 
both cases, making hi^q 2 ^) p» h na i ve (q\). For smaller q±, on the other hand, the naive 
expectation gives a significant overestimate. This can also be concluded from Fig. pi], 
which shows the distribution in transverse energy, dE '±/ 1 dydq\ = q± ■ dn/dydq\. We see 
that the total E± in the minijet region, q± ~ 5 GeV, is almost a factor 2 larger in the 
naive estimate. It is often realized that this large contribution for very low qj_ must be 
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0.5 1 5 10 50 100 500 

<?JGeV) 

Figure 9: The ratio h(qj_) / '/2. n cmie(<?i_) between the cross section in our approach and the 
naive expectation in eq. §2t\) , as a function of q±. (Note the logarithmic scale on the 
x-axis.) 



unphysical, and therefore a soft cutoff, q±, m i n , is introduced in many calculations. This 
cutoff is often assumed to be around 2 GeV and slowly growing with energy (see e.g. 0). 
In Fig O we show the integrated transverse energy f Qj - dq\ ^ E ,\ . From this figure we see 



that 



a dE - 



Ql dydq 



, o dE 
dq 

ourresult " y X,min 



- ,J ~ dydqj 



(29) 



for q±.min ~ 2.1 GeV. (Remember that the low fcj_-cutoff Qq was given by 0.6 GeV.) This 
means that for asymptotic energies the total E± in our approach equals the naive result 
when q±, m in ~ 2.1 GeV. For smaller energies the maximal value of q± is limited, which 
implies that a somewhat smaller value of q±, m in is needed. A reliable quantitative estimate 
would need a calculation which includes quark jets and non-leading contributions, which 
could be obtained with the help of the MC program. We note in particular that in our 
approach the corresponding effective cutoff of the naive approach, q± jm i n , saturates for 
very large energies. In the conventional approach it is difficult to make predictions for 
higher energies without a physical understanding of the energy dependence of qx,min- 

Only for comparison, we have in Fig. |Tt] also included the distribution one would 
obtain from eq. ( p5| ) if the structure function F were scaling and independent of q]_. (The 
normalization is adjusted to our result for large momenta.) As is seen, this would give a 
very much larger (and totally unrealistic) increase for small q±. 

We end this section by noting that one chain forms a set of correlated jets. The fact 
that the /^-distribution saturates for long chains (cf eq. flU])) implies that not only the 
p± of the jets, but also their density (the number of jets per unit rapidity) is independent 
of the length of the chain, i. e. independent of the energy in the collision. The jet 
density grows because the number of possible chains increases proportional to s x . It is 
also conceivable that the number of chains is not random, but e. g. correlated with 
the impact parameter, so that central collisions have more and peripheral collisions fewer 
chains [10[|. In this paper we have only studied the average jet multiplicity, and we 



postpone the study of different types of correlations to future investigations. 
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Figure 10: Transverse energy distribution, dE±/dydqj_, according to our calculations 
(continuous curve) and the "naive" estimate in eq. ([2%) (dashed curve). For comparison 
we also show the result from scaling structure functions, F(x) independent of q\ (dotted 
curve). The normalization is arbitrary. 

The possible correlation between different chains may also affect the total cross section. 
We want to stress that the results in eq. fl25|) correspond to the cross section for jet 
production. To find the number of jets per event we have to divide by the total cross 
section, which is not directly obtained from our formulae. 

6 Conclusions 

A good understanding of jet and minijet production, and transverse energy flow, is essen- 
tial for a proper interpretation of new phenomena in pp collisions at LHC or a possible 
formation of a quark-gluon plasma in nucleus collisions at RHIC or LHC. For very high pj_ 
the parton flux can be described by DGLAP (/^-ordered) evolution. The large momentum 
transfer between the colliding partons implies that the two evolving chains are indepen- 
dent of each other. For moderate and smaller p± (and high energies) we enter the BFKL 
regime and non-ordered chains become important. Here it is essential to take coherence 
effects and correlations into account. As we can then have several hard subcollisions in a 
single chain or ladder diagram, it is also important to avoid double counting. 

The parton evolution in the BFKL regime can be described e. g. by the CCFM model 
or the semisoft formalism. The LDC model is a reformulation and generalization of the 
CCFM formalism, and the symmetric structure of the LDC model makes it particularly 
suited for a description of (mini)jet production. The result can be interpreted as the 
production probability for exclusive final states, which makes it convenient for treating 
the problem of double counting. The (mini)jet cross section is described in a /^-factorizing 
form in terms of non-integrated structure functions and off-shell subcollision cross sections. 
The result is a dynamical suppression of small-p_i_ jets, which removes the strong sensitivity 
to a low-p_i_ cutoff in "naive" estimates based on integrated structure functions. 
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Figure 11: Integrated transverse energy distribution, J 9± dq±dE±/dydq±, according to 
our calculations (continuous curve) and the "naive" estimate in eq. j[2t\) (dashed curve). 
The normalization is as in Fig. [77]. 
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8 Appendix 

In this appendix we derive the explicit form of the function h(qj_), which appears in 
the expression for the total cross section in eq. fl25|) , for the case of a running coupling 
constant a s oc 1/ ln(Ar[/A 2 ) = 1/ k. In each vertex the scale in a s is taken to be the largest 
transverse momentum, which due to the consistency constraint in eq. @ approximately 
coincides with the largest virtuality. Thus we write the cross section given by eq. ([H]) in 
the form 

g§^-. 2 tei)-%i>- (30) 

As before, we consider separately each of the three possible cases for the relative sizes of 
the transverse momenta of the gluons k a , k, and kb in Fig. [| 

Case 1. k± > k±^ a , k± t b, which implies that q± yCL ~ q± t b ~ k±, 

Case 2. k± t b > k± > k± >a , which implies that q± a « k± and qj_,b ~ k± jb , 

Case 3. k± < k± >a , kj_,bi which implies that q± A ~ k^ a and q± t b ~ k± ; b- 

We also assume the power-like approximation for f(n) in eqs. ( |I^|T5| ), T{x, k) ~ 
x~ x ■ f(n) oc x~ x k^ 1 ~ 1 ~ x~ x ■ k 2 , where ^ — 1 is approximately equal to 2. 

We start by considering the contribution from case 1. Here the scale in a s is given 
just by k\ q\. Thus the only modification to eq. (pi]) is a factor a 2 s {q\) oc 1/ ln 2 ^), 



which is factored out explicitly in the definition of h(qj_) in eq. (30). Thus we find 



da 
dq\ 



oc 2 



contr 1 



1 

S 'h? 



(31) 



We have a factor 2 in the inclusive cross section, because we count both outgoing 
partons. When this expression is inserted into eq. (pOD, the integrations with respect to 
k^ a and k L)h give 



dn 



«0 



K 



1 

P 



K„ ■ Kl 



f 2 



1 

9 



3\2 



(32) 



where l q = lng^ and k = IuQq. Here the factor 1/2 compensates for double counting 
as described in the main text. The remaining integrations with respect to x a and Xb (cf 
eq. (|23"D) give finally: 



da 



dqj_dy 



•oc 



contr 1 



<zl +2A 



1. 1. 

Z 2 A 9 9 



3N2 



(33) 
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For case 2, the subcollision cross section ^p- is given by 



da 



dq\ 



OC 



contr 2 



J^ b dkl ■ [5(ql - kl b ) + 5(ql - kl)} ■ J- • 1 



rk 2 




j K ± 


'"dk 2 ± - 


Jk 2 


1 


1 


k 2 


K b 



1 1 



%! - In(-) + 0(kV b - - tfj • -4- • 7- (3 i) 



After insertion into eq. (pD[) , calculations analogous to the ones in case 1 give the 
following result for the contribution to the total cross section: 



da 



dq]_dy 



oc 



contr 2 



1 1 



qi +2X l\ 3 
+^g ' (lq ~~ K 0y 



A 



76 K 3/3 

3 ~ ~3~ 

l + la 



V 



+ 



A(A + 1) (A + l) : 



In case 3 the appropriate expression for the subcollision cross section Jp- is 



(35) 



da 



dq\ 



oc 



case 3 



[ kla dkl ■ [6(ql - klj + 5(ql - k%)\ ■ -J- • J- • 1 • 1, (36) 



where the lower integration limit is given by max 
and performing the integrations we obtain 



k 2 k 2 



6 Ql). Inserting this into eq. (|20|) 



<icr 



dqldy 



oc 



contr 3 



<zi +2A 



A(A + 1) 



- «o \ 

■ ^—Z h L - «0 ~ 



(A + l) ; 



• ^g - (^g — ^0) 



(37) 



Adding the three contributions, and including the symmetric situation k± ta > &±,& to 
case 2, we find finally for the function h(ql) in eq. (|30|) 



+ 



+ 



A 9 v 






1 1 




3/3 


A3 
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3 


1 




Z 2 

'« l 2 


A(A + 1) 



2 



V 

K , 



+ g ' ' (^g ~~ K o) 



l + l Q 



A(A+1) (A + l) 1 



(a + i; 



2 1 



la ■ (la -Ko)>- 



This can be compared with the "naive" expression obtained from eq. With the same 
normalization we get 



1 / I 2 



(A + l) 2 (A + l) 



(39) 



where the factor 2 in the inclusive cross section follows because each hard subcollision 
produces two jets. We see that the leading term for large q±, ^(lng^) 6 , is the same in 
both cases, but h na i Ve is significantly larger for smaller g^-values. 
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